CoefRcients of bosonized dimer operators in spin-| XXZ chains and their applications 
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we determine coefficients d^y and dz of 

-iy(s]S]+^ + 

We also calcu- 



Comparing numerically evaluated excitation gaps of dimerized spin 
formula for the low-energy effective sine- Gordon theory, 

bosonized dimerization operators in spin-i XXZ chains, which are defined as 
SlSj+i) = d^ySm{V^(t>{x)) + --- and (-iyS'|S|+i = d,,sin{.j^4>{x)) + ■ ■ ■ 
late the coefficients of both spin and dimer operators for the spin-i Heisenberg antiferromagnetic 
chain with a nearest-neighbor coupling J and a next-nearest-neighbor coupling J2 = 0.2411J. As 
applications of these coefficients, we present ground-state phase diagrams of dimerized spin chains 
in a magnetic field and antiferromagnetic spin ladders with a four-spin interaction. The optical 
conductivity and electric polarization of one-dimensional Mott insulators with Peierls instability are 
also evaluated quantitatively. 
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I. INTRODUCTION 

Quantum magnets in one dimension are a basic class 
of many-body systems in condensed matter and statisti- 
cal physics (see e.g., Refs. |l|,||). They have offered various 
kinds of topics in both experimental and theoretical stud- 
ies for a long time. In particular, the spin-i XXZ chain is 
a simple though realistic system in this field. The Hamil- 
tonian is defined by 

ti = J y (j,- 'J-,-1-1 + Ojb^j^i -f- AzSjSjj^i), (1) 



where 5" is a-component of a spin- 2 operator on j- 
th site, J > is the exchange coupling constant, and 
A^ is the anisotropy parameter, ^h^^ model is ex- 
actly solved by integrability methodsj3a and the ground- 
state phase diagram has been completed. Three phases 
appear depending on A^; the antiferromagnetic (AF) 
phase with a Neel order (S"!) = -(5*1+1) (A^ > 1), the 
critical Tomonaga-Luttinger liquid (TLL) phase (—1 < 
A^ < 1), and the fully polarized phase with {S^-) = 1/2 
(A2 < — 1). In and around the TLL phase, the low- 
energy and long-distance properties can be understood 
via effective field theory techniques-sjich,as bosonization 
and conformal field theory (CFT)JdMBtl These theoret- 
ical results nicely explain experiments of several quasi 
one-dimensional (ID) magnets. The deep knowledge 
of this model is also useful for analyzing plentiful re- 
lated magnetic systems, such as spin-g chains with some 
perturbations (e.ft_ j— .external fields H additional mag- 
netic anisq tropies,Erllj dimerizationlljllj), coupled spin 
chainsiSEjZl spatially anisotropic 2D or 3D spin sys- 
tems Jl3c3 etc. 

A recent direction of studying spin chains is to es- 
tablish solid correspondences between the model (|l|) and 
its effecti\je-|i|h*ory. For example, Lukyanov and his col- 
laboratorg22r£3 have analytically predicted coefficients of 



bosonized sp>i«|-aperators in the TLL phase. Hikihara 
and Furusakicjl23 have also determined them numerically 
in the same chains with and without a uniform Zeeman 
term. Using these results, one can now calculate ampli- 
tudes of spin correlation functions as well as their critical 
exponents. Furthermore, effects of perturbations on an 
XXZ chain can also be calculated with high accuracy. 
It therefore becomes possible to quantitatively compare 
theoretical and experimental results in quasi ID mag- 
nets. The purpose of the present study is to attach a 
new relationship between the spin-i XXZ chain and its 
bosonized effective theory. Namely, we numerically eval- 
uate coefficients of bosonized dimer operators in the TLL 
phase of the XXZ chain. Dimer operators {—lyS'^S'^j^i, 
as well as spin operators, are fundamental degrees of free- 
dom in spin- 2 AF chains. In fact, the leading terms of 
both bosonized spin and dimer operators have the same 
scaling dimension 1/2 at the S'C/(2)-symmetric AF point 
A^ = 1 (see Sec. ||). 

In Rcfs. |2^ , p5| , Hikihara and Furusaki have used 
density-matrix renormalization-group (DMRG) method 
in an efficient manner in order to accurately evaluate co- 
efficients of spin operators of an XXZ chain in a magnetic 
field. Instead of such a direct powerful method, we utilize 
the relationship between a dimejrized XXZ chain and its 
effective sine- Gordon theorylldcj to determine the coeffi- 
cients of dimer operators (defined in Sec. ||), i.e., excita- 
tion gaps in dimerized spin chains are evaluated by nu- 
merical diagonalization method and are compared with 
the gap formula of the effective sine- Gordon theory. In 
other words, we derive the information on uniform spin- 
i XXZ chains from dimerized (deformed) chains. More- 
over, we also determine the coefficients of both spin and 
dimer operators for the spin-i Heisenberg (i.e., XXX) AF 
chain with an additional next-nearest-neighbor (NNN) 
couplin g J2 — 0.2411J in the similar strategy. As seen 
in Sec. HID, evaluated coefficients are more reliable for 
the J-J2 model, since the marginal terms vanish in its 



effective tlieory. 

The plan of this paper is as follows. In Sec. ||, we 
shortly summarize the bosonization of XXZ spin chains. 
Both the XXZ chain with dimerization and the chain in 
a staggered magnetic field are mapped to a sine-Gordon 
model. We also consider the AF Heisenberg chain with 
NNN couphng J2 = 0.2411 J. In Sec. ffl, we explain how 



to obtain the coefficients of dimer and spin operators by 
using numerical diagonalization method. The evaluated 
coefficients are listed in Tables | and || and Fig. |[ These 
are the main results of this paper. For comparison, the 
same dimer coefficients are also calculated by using the 
formula of the ground-state energy of the sine-Gordon 
model. We find that the coefficients fixed by the gap 
formula are more reliable. We apply these coefficients to 
several systems and physical quantities related to an XXZ 
chain (dimerized spin chains under a magnetic field, spin 
ladders with a four-spin exchange and optical response 
of dimerized ID Mott insulators) in Sec. |V[ Finally our 
results are summarized in Sec. M. 



II. DIMERIZED CHAIN AND SINE-GORDON 
MODEL 



can be exactly evaluated from Bethe ansatz 
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Here we have introduced new parameters 77 and R. The 
former is the critical exponent of two-point spin corre- 
lation functions and used in the discussion below. The 
latter is called the compactification radius. It fixes the 
periodicity of fields 4> and 6 as (J)/Vk ~ (I)/VK + 27ri? 
and ^fKO ~ ^fK9 -t- \/R. Using the scalar fields (j) and 0, 
we can obtain the bosonized representation of spin oper- 
ators: 
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dx(t) + (-l)^ai cos(V47r(?!)) -|- 
boi-iy + bi cos(V4^(/)) 4- 



(5a) 
(5b) 



where a„ and 6„ are non-universal constants, and some 
of them with small n have been determined accurately 
in Refs. pl^-pa_. la this formalism, vertex operators are 
normalized asEircj 



In this section, we explain the relationship between a 
dimerized XXZ chain and the corresponding sine-Gordon 
theory in the easy-plane region — 1 < A^ < 1. XXZ 
chains in a staggered field and the AF Heisenberg chain 
with NNN coupling J2 — 0.2411J are also discussed. The 
coefficients of dimer operators are defined in Eq. (R). 



A. Bosonization of spin- 5 XXZ chain 

We first review the effective theory for undimerized 
spin chain (|]). According to the standard strategy, XXZ 
Hamiltonian (nl) is bosonized as 



nfi^ = / dx{^[K-\dx<pf+Kidx9f 



/Qi<]<l>{x)^-iq4>{x') 



at |a;- a;'| > a. (6) 



This means that the operator e'^"**^^ 
sion Kq'^/{4:n). 



has scaling dimen- 



In addition to the spin operators, the-basionized forms 



of the dimer operators are known to b 



s 



(2) 



i-iy{S^S^+, + Sp^^,) ^dxy sin(\/4;^0) + • • • , (7a) 
(-l)^S'|S'|+i «4 sin(V4^0) + ■■■ . (7b) 

In contrast to the spin operators, the coefficients dxy and 
dz have never been evaluated so far. To determine them 
is the subject of this paper. It seems to be possible to 
calculate dx^j.z by utilizing Ii'jln(|5|) and operator-product- 
expansion (OPE) technique ,ErLl buL.it requires the correct 
values of all the factors a„ and 6n.L3 Therefore, we should 
interpret that the dimer coefficients dxy,z are independent 
of spin coefficients a„ and 6„. 



in the TLL phase. Here, (j){x) and 9{x) are dual scalar 
fields, which satisfy the commutation relation. 



[cj){x),9ix')] = -iAtcpix~x'), 



(3) 



with X — ja (a is the lattice spacing). As we see in 
Eq. (||), cos{VT6tt4>) is irrelevant in — 1 < Az < 1, and 
becomes marginal at the S'C/(2)-symmetric AF Heisen- 
berg point A^ = |l. || T he coupling constant A has been de- 
termined exactlyB30 Two quantities K and v denote the 
TLL parameter and spinon velocity, respectively, which 



B. Bosonization of dimerized spin chain 

Next, let us consider a bond- alternating XXZ chain 
whose Hamiltonian is given as 

^xxz-s ^jj2 [(1 + {-iys,y)iSJSJ+, + spl,) 



+(A, -f i-iys^jsp^^,] 



(8) 



In the weak dimerization regime of |<5a;y^2| <C 1, the 
bosonization is applicable and the dimerization terms can 



be treated perturbatively. From the formula (0), the ef- 
fective Hamiltonian of Eq. (| 
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o/XXZ-^ 
^eff 



dx 



+ -{d^yd^y + S,d,)sm{V4^(j))^ |. (9) 

Here, we have neglected all of the irrelevant terms in- 
cluding cos(vT67r(/)). This is nothing but an integrable 
sine-Gordon model (see e.g., Refs. pl]j2q and references 
therein). The sin(-\/47r(/)) term has a scaling dimension 
K, and is relevant when K < 2, i.e., -0.7071 < A^ < 1. 
In this case, an excitation gap opens and a dimerization 
{S^Sf^^ - S^^-^^Sf^^) "^ occurs. The excitation^pec- 
trum of the sine-Gordon model has been known,OE3 and 
three types of elementary particles appear; a soliton, the 
corresponding antisoliton, and bound states of the soli- 
ton and the antisoliton (called breathers). The soliton 
and antisoliton have the same mass gap Eg. There ex- 
ist [477 — 1] breathers, in which \A\ stands for the integer 
part of A. The mass of soliton and n-th breather Eb^ 
are related as follows. 



Eb-^ = 2Es sin 



2(47?-!) 



n — 1, 



, [4^ - 1] 



(10) 

The breather mass in units of the soliton mass is shown 
in Fig. n^ as a function of A^. Note that there is no 
breather in the ferromagnetic side Az < 0, and the 
lightest breather with mass Eb^ is always heavier than 
the soliton in the present easy-plane regime. Following 
Refs. 2llp9|, the soliton mass is also analytically repre- 
sented as 
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(11) 



In addition, the difference between the ground-state en- 
ergy ffi-ee of the free-boson theory (g) with A = per site 
and that of r-Ujp, sine- Gordon theory (|^), £sGj has been 
predicted asEiJE3 
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However, we should note that the above formula is invalid 
for the ferromagnetic side A^ < (77 < 1/2) since it 
diverges at the XY point A^ = (77 = 1/2). 

A similar sine-Gordon model also emerges in spin-^ 
XXZ chains in a staggered field. 



^stag ^ ^XXZ ^ ^(_1)J7j^5'| 



(13) 




FIG. 1: (Color online) Ratio of the 7i-th breather mass _Es„ 
to the soliton mass Es as a function of the XXZ anisotropy 
Az in the sine-Gordon model (13) or wM- 



The staggered field hs induces a relevant perturbation 
cos{-\f^(j)). Therefore, the resultant effective Hamilto- 
nian is 



.T/stag -T/XXZ 



/ Ax—ai cos( 
J a 



47r0). (14) 



If we redefine the scalar field 4> as (j) + V^/ii the form of 
Eq. (|l4| ) becomes equivalent to that of Eq. (g). Thus, the 
soliton gap of the model (Oh is equal to Eq. (O) with 
the replacement of Sxyd^y + Szdz — >■ hgai/J. Namely the 
soliton gap of the model (O) is given by 
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(15) 



This type of staggered-field induced gaps has been ob- 
served in some quasi ID magnets with an alternating gy- 
romagnetic tensor oij-|PsAi^loshinskii-Moriya interaction 
such as Cu benzoate.crt3e3 

Masses of the soliton, antisoliton and breathers are re- 
lated to the excitation gaps of the original lattice sys- 
tems, Eqs. (||) and (|3|). The soliton and antisoliton 
correspond to the lowest excitations which change the 
z component of total spin S^^^ = J2-j ^j by ±1- On the 
other hand, the lightest breather is regarded as the low- 
est excitation with AS'^^^ = 0. At the 5[/(2)-symmetric 
AF point A^ = 1, there are three breathers. The soli- 
ton, antisoliton and lightest breather are degenerate and 
form the spin-1 triplet excitations (so-called magnons). 
The second lightest breather is interpreted as the sin- 
glet excitation with A^tot = 0. In the ferromagnetic 
regime A^ < 0, where any breather disappears, the low- 
est soliton-antisoliton scattering state would correspond 
to the excitation gap in the sector of AS'f^^ = 0. 



C. J- J 2 antiferromagnetic spin chain 

In the previous two subsections, we have completely 
neglected effects of irrelevant perturbations in the low- 
energy effective theory. However, as already noted, the 
A term becomes nearly marginal when the anisotropy A^ 
approaches unity. In this case, the A term is expected to 
affect several physical quantities. Actually, such effects 
have been studied in both the models (j^) [Ref. |l^ and 
([|) [Refs. |M. 

It is knownl3 that a small AF NNN coupling J2 de- 
creases the value of A in the S'C/(2)-jSi'mmetric AF Heisen- 
berg chain. Okamoto and NomuraEil have shown that the 
marginal interaction vanishes, i.e., A — ^ in the following 
model: 



H" 



^^{JSj ■ Sj+i + J2SJ ■ s 



J+2) 



(16) 



with J2 = 0.241 IJ. On the J2/J axis, this model is lo- 
cated at the Kosterlitz-Thouless transition point between 
the TLL and a spontaneously dimerized phase. From 
this fact, if we replace H^^^ with H""" in the SU{2)- 
symmetric models (pi) and (13), namely, if we consider 
the following models: 
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XXX-(5 
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^^nnn ^ ■^(_ 1)J,5 J^^. . S^ + ^, (17a) 
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=H""" + ^(-l)^/is^|, (17b) 



n 



then their effective theories are much closer to a pure 
sine- Gordon model. In other words, the predictions from 
the sine-Gordon model, such as Eqs. (11) and (|l5|), be- 
come more reliable. 



III. COEFFICIENTS OF DIMER AND SPIN 
OPERATORS 

From the discussions in Sec. |y, one can readily find a 
way of extracting the values of dxy,z and ai in Eqs. (M) 
and (|^) as follows. We first calculate some low-energy 
levels in S^^^. = ±1 and S^^^. = sectors of the models 
(§); (JH) and (|7|) by means of numerical diagonaliza- 
tion method. Since all the Hamiltonians (pi), (Ol) and 
( pTI ) commute with S^^^ — J2j ^j ' the numerical diago- 
nalization can be performed in the Hilbert subspace with 
each fixed S'f^^. In order to extrapolate gaps to the ther- 
modynamic limit with reasonable ajCcucacy, we use ap- 
propriate finite-size scaling methods&Ej for spin chains 
under periodic boundary condition (total number of sites 
L = 8, 10, • • • , 28, 30). Secondly, the coefficients dxy,z 
and ai of the spin-^ XXZ chain and the J-J2 chain are 
determined via the comparison between the sine- Gordon 
gap formula (|Tl|) and numerically evaluated spin gaps for 
various values of Sxy,z and hg. In this procedure, (as al- 
ready mentioned) the energy difference between the low- 
est (i.e., ground-state) and the second lowest levels of 



the S'^ot = sector (gap with AS^^^ = 0) and that be- 
tween the ground-state level and the lowest level of the 
5fot = ±1 sector (gap with AS^^^ — ±1) are respec- 
tively interpreted as the breather (or soliton-antisoliton 
scattering state) and soliton masses in the sine-Gordon 
scheme. 



A. TLL phase and Numerical diagonahzation 

In this subsection, wc focus on the TLL phase of uni- 
form spin-i XXZ chains (|l|) and test the reliability of our 
numerical diagonahzation. The low-energy properties are 
described by Eq. (g), which is a free boson theory (i.e., 
GET with central charge c = 1) with some irrelevant per- 
turbations. Generally, the finite-size scaling formul a for , 
the excitation spectrum in any GET has been provedtHlLJ 
to be 



AEc 



Ec 



27n; , , 



(18) 



Here Eq and Eq are respectively the ground-state energy 
and the energy of an excited state generating from a pri- 
mary field O in the given GET. Remaining quantities 
[O], V, and La are the scaling dimension of the operator 
O, the excitation velocity and the system length, respec- 
tively. In the case of the spin chain (|l|), the bosonization 
formula (|^) indicates that E^±i^e and E^±i^/j:^^ corre- 
spond to the excitation energies in the S^^^ = ±1 and 



qz 



sectors, respectively. The irrelevant pertur- 



bations can also contribute to the finite-size correction 
to excitation energies. From the C/(l) and translational 
symmetries of the XXZ chain (p, one can show that 
the finite-size gap AE/^s" =±1 has no significant modi- 
fication from the perturbations, while the correction to 

AE^s^^t=Q is proportional to L^~[°' " 1. Therefore, the 
following finite-size scaling formulas are predicted: 

27rt; 1 

AEas' -±i ~ \ , 

^itot-±i La AK 
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(19a) 
(19b) 



with Co being a non-universal constant. Here we have 



used [e 



ln^/^T9^ 



n^l(AK) and [e 



in\/47rd 



n^K. At 



the 5't/(2)-symmetric AF point A^ = 1, A^ASt-„t=±i = 
AE^S' =0 = A£'su2 holds and the marginal A term mod- 
ifies the scaling form of the spin gap. The marginal texm 
is known to yield a logarithmic correction as follows:u 



AP ^ 27rw / 1 ci C2 

^^^"' ^ L^ U 1^ (h^ 



(20) 



Here ci,2 are non-universal constants. 

As an example, numerically evaluated gaps with 
AS't^^t = ±1 and AS't^ot = in the case of Az = 0.6 
are respectively represented as circles and triangles in 
Fig. g(a). Gircles are nicely fitted by the solid curve 
AEas;,=±i/J = 8.019 x 10"'' + 2.977/L. This result 



(a) 0.7 

^0.6 


§,0.5 


|0.4 
|0.3 

:lo.2 

0.1 






'. A2=0.6 




A 






A 






a" ^ 


~ 


A^' 


^^'. 


: x^ 


r^ 


^^^^^ ; 






A aEas^„=o/-.' ■ 


'^^^^^""^ . 




• aEas^„,= ±i/J- 



0.05 0.1 

1 //. 



0.15 



(b) 



Q.0 

cc 

S'o 

u5 

I 

■fO, 



5 

A 


: Az= 


' ' 1 
=1 




' 1 ' 


X 


^ , 


3 


- 




X' 






- 


2 


- 


y^ 


r 






- 


1 


- . 


/^ 


• 


AEsu2 


IJ 


- 







. . 1 




. 1 . 




- 



0.05 0.1 

1 //. 



0.15 



FIG. 2: (Color online) (a) Numerically evaluated gaps with 
AS'iot = ±1 (circles) and ASf^t = (triangles) for XXZ 
chains with A^ = 0.6 and finite length L. The solid curve 
8.019xl0"'* + 2.977/L (dashed curve 1.312 x 10"^ + 5. 982/L- 
4.764/L^'*^'^®) is determined by fitting the circles (trian- 
gles), (b) Gaps of finite-size Heisenberg chains with A^ — 1. 
The solid curve is AE,u2/J = 2.173 x 10"" + 4.965/L - 
2.203/(LlnL) + 1.200/(L(lnL)2). 



is consistent with the fact that an easy-plane anisotropic 
XXZ model is gapless in the thermodynamic limit and 
that the exact coefficient of the 1/L term is 2ttv/{4:JK) = 
3 at Az = 0.6. Similarly, triangles can be fitted by 



AK 



ASt-ot^ 



=o/J 



1.312 X 10-3 + 5.982/L - 4.764/Li-8376 
where 1.8376 = 1- iK. The factor 5.982 of the 1/L 
term is very close to 2'KvK/{Ja) = 6.040. The spin gap 
at 5'C/(2)-symmetric point is also represented in Fig. g(b). 
Following the formula (^0|), we can correctly determine 
the fitting curve AEsu2/J = 2.173 x lO"" + AMh/L - 
2.203/(LlnL) -I- 1.200/(i(lnL)2), in which the factor of 
the second term is nearly equal to Trv/{Ja) = 4.935. 
These results support the reliability of our numerical di- 
agonalization. We note that a more precise finite-size 
scaling analysis for AF Heisenberg model has been per- 
formed in Ref. Udi 



B. Dimer coefHcients of XY model 
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FIG. 3: (Color online) Numerically evaluated gaps 
A-Eas^ =±i (circles) and A_Bas' =o (triangles) in dimerized 
XY models with A^ = Jz = 0. Solid and dashed lines are 
the exact results determined via Jordan- Wigner transforma- 
tion. These lines respectively correspond to the soliton and 
breather masses in the framework of sine-Gordon theory. 



In this subsection, we consider a special dimerized XY 
chain with A^ = Sz = 0. It is mapped to a solvable 
free fermion system through Jordan- Wigner transforma- 
tion. Therefore, our numerically determined coefficients 
in Eq. (|7|) can be compared with the exact value. The 
lowest energy gap with AS^^^ = ±1, which corresponds 
to the soliton mass Es, is exactly evaluated as 



A.EASf„,=±i/^ 



m) 



Comparing Eq. ( |2l| ) with Eq. (^jj), we obtain the exact 
coefficient 



d^y = I/tt = 0.3183 



(22) 



at the XY case A, 



0. The exact solution also 



tells us that the excitation gap with AS'^^ 



is 
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^o/J = '2Sx 



(23) 



This is consistent with the sine- Gordon prediction that 
any breather disappears and the relation Ebi = ^.Es 
holds just at the XY point Az = 0. 

Figure |^ shows the comparison between the energy 
gap calculated by numerical diagonalization with Aitken- 
Shanks process and Eq. ( |2l] ) [or Eq. (p3|)]. Except for 
AE^S' =0 in the weak dimerized regime S^y ^ 0.1, nu- 
merically calculated gaps coincide well with the exact 
value. We have found that when Sxy^z becomes smaller, 
the precision of Aitken-Shanks method is decreased due 
to a large size dependence of gaps. 



Next, let us move onto the evaluation of excitation 
gaps in dimerized XXZ chains. In this case, since the 
system is not critical, the above finite-size scaling based 
on CFT cannot-be applied. Instead, we utilize Aitken- 
Shanks methodcS to extrapolate our numerical data to 
the values in the thermodynamic limit. 



C. Dimer coefficients of XXZ model 

In the easy-plane region — 1 < A^ < 1, any generic 
analytical way of determining the coefficients in Eq. (1^) 
has never been known except for the above special point 
Az = Sz = 0- To obtain dxy (respectively dz), we 



TABLE I: Dimer coefficients {dxy and dz), TLL parameter K, compactification radius R, spinon velocity v of spin- j XXZ 
cliain. Dimerization-induced gaps are also listed in the final column. The final line is the result for the J-J2 chain (|l6|). The 
same data of dxy,z are also shown in Fig. H. 
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0.9 
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J-J2 model 



0.228 
0.278 
0.297 
0.309 
0.318 
0.324 
0.327 
0.328 
0.328 
0.324 
0.318 
0.309 
0.297 
0.278 
0.252 
0.213 
0.364 



(0.204) 
(0.261) 
(0.284) 
(0.299) 
(0.310) 
(0.318) 
(0.323) 
(0.325) 
(0.325) 
(0.323) 
(0.318) 
(0.311) 
(0.302) 
(0.289) 
(0.273) 
(0.248) 
(0.361) 



0.110 
0.141 
0.154 
0.165 
0.174 
0.182 
0.188 
0.193 
0.197 
0.200 
0.202 
0.202 
0.200 
0.194 
0.184 
0.163 
0.188 



(0.097) 
(0.131) 
(0.146) 
(0.159) 
(0.169) 
(0.177) 
(0.185) 
(0.191) 
(0.196) 
(0.200) 
(0.203) 
(0.204) 
(0.204) 
(0.203) 
(0.199) 
(0.191) 
(0.182) 



0.5 

0.5838 

0.6288 

0.6695 

0.7094 

0.75 

0.7924 

0.8375 

0.8864 

0.9401 
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1.068 

1.147 

1.241 

1.355 

1.5 

0.5 



0.3989(= 1/V2^) 

0.3692 

0.3557 

0.3448 

0.3349 

0.3257 

0.3169 

0.3082 

0.2996 

0.2910 

0.2821(= 1/V^) 

0.2730 

0.2634 

0.2533 

0.2423 

0.2303 

0.3989 



1.571(= n/2) 

1.518 

1.465 

1.410 

1.355 

1.299 

1.242 

1.184 

1.124 

1.063 

1 

0.9353 

0.8685 

0.7990 

0.7263 

0.6495 

1.174 
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0.6667 

0.7061 
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0.8281 
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1.550 

2 

0.6667 



numerically calculate excitation gaps at the points 5xy 

Al- 



0.35 



(<5. 



0.05, 0.1, •••, 0.3 with fixing 5^{5xy) = 0. 



though both Ai?ASf„t=o and Ai?ASf„t=±i are applicable 
to determine d^y^z hi principle, we use only the latter gap 
since it more smoothly converges to its thermodynamic- 
limit value via Aitken-Shanks process, compared to the 
former. In fact, Eq. (19) suggests that AE^s^ =o is 



subject to effects of irrelevant perturbations and there- 
fore contains complicated finite-size corrections. Coeffi- 
cients dxj, (dz) can be determined for each ^™ (Sz) from 
Eq. (|ll|). Since the field theory result (^ is gener- 
ally more reliable as the perturbation S^y^z is smaller, 
we should compare Eq. (11) with excitation gaps deter- 
mined at sufficiently small values of Sxy,z- However, the 
extrapolation to thermodynamic limit by Aitken-Shanks 
method is less precise in such a small dimerization region 
mainly due to large finite-size effectsJlil'Lj Therefore, we 
adopt coefficients dxy,z extracted from the gaps at rela- 
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FIG. 4; (Color online) XXZ-anisotropy (A^) dependence of 
dimer coefficients dxy and dz. Filled [Open] circles repre- 



tively large dimerization 4y(^) = 0.1 and 0.3, and they gent d^y determined from dimerization gap at (<5^a,(5; 
are listed in Table |: the values outside [inside] parenthe- 



ses are the data for Sxy(z) = 0.3 [0.1]. The anisotropy 
dependence of the same data dxv,z is depicted in Fig. ^. 
The data in Table | and Fig. H are the main result of 
this paper. The difference between dxy{z) outside and in- 
side the parentheses in Table | could be interpreted as the 
" strength" of irrelevant perturbations neglected in the ef- 
fective sine-Gordon theory or the "error" of our numerical 
strategy. The neglected operators must bring a renor- 
malization of coefficients dxy,z, and the "error" would 
become larger as the system approaches the Heisenberg 
point since (as already mentioned) the A term becomes 



(0.3,0) [— (0.1,0)]. Similarly, filled [open] triangles show 
dz determined from dimerization gap at {Sxy,Sz) = (0,0.3) 
[=(0,0.1)]. 



marginal at the point. 

We here discuss the validity of the numerically deter- 
mined dxy,z in Table | and Fig. ^. Table | shows that in 
the wide range —0.3 < A^ < 0.9, the difference (error) 
between dxy.z outside and inside the parentheses is less 
than 8 %. As expected, one finds that the error gradu- 
ally increases when the anisotropy Az approaches unity. 
Similarly, the error is large in the deeply ferromagnetic 



regime A^ < —0.3. This is naturally understood from 
the fact that as A^ is negatively increased, the dimer- 
ization term sin(\/47r0) becomes less relevant and effects 
of other irrelevant terms is relatively strong. Indeed, for 
Az < -0.7071 (K > 2), the dimerization does not yield 
any spin gap and our method of determining d^y^z cannot 
be used. Furthermore, it is worth noting that the spin 
gap is convex downward as a function of dimerization 
5xy,z in the ferromagnetic side A^ < 0, and the accuracy 
of the fitting therefore depreciates. 

In addition to coefficients d^y^z, let us examine dimer- 



(a) 0.7. 



ization gaps and the quality of fitting by Eq. ( |ll| ) 
citation gaps for A^ = 0.6 are shown in Fig. 
an example. Remarkably, both soliton-gap curves 
with the values d^y^z outside and inside the parenthe- 
ses in Table I fit the numerical data AEas^^^=±i in the 
broad region < Srcy{z) < 0.3 with reasonable accuracy. 
The former solid curve is slightly better that the latter. 
The breather gaps AE^s' =o a-nd corresponding fitting 
curves are also shown in Fig. H. This breather curve is 
determined by combining the solid curve ( ]ll| ) and the 
soliton-breather relation dl^). It slightly deviates from 
numerical data, especially, in a relatively large dimeriza- 
tion regime 0.15 < Sxy(z)- As mentioned above, this de- 
viation would be attributed to irrelevant perturbations. 
The breather-soliton mass ratio Eb^ /Eg [see Eq. (|l(]|)] in 
the sine-Gordon model (0) and the numerically evaluated 
Ai?AS^ =o/ AE/^S' =±1 are shown in Fig. g. These two 
values are in good agreement with each other in the wide 
parameter region —0.5 < A^ < 1, although their differ- 
ence becomes slightly larger in the region 0.5 ^A^ < 1, 
which includes the point A^ = 0.6 in Fig. g. Gaps 
Ai?AS^ =±1 for dimerized XXZ chains with several val- 
ues of both 5j:y and 5z are plotted in Fig. 0. It shows 
that the numerical data are quantitatively fitted by the 
single gap formula (Jll|). All of the results in Figs. |^0 in- 
dicates that a simple sine-Gordon model ^ can describe 
the low-energy physics of the dimerized spin chain (0) 
with reasonable accuracy in the wide easy-plane regime. 
This also supports the validity of our numerical approach 
for fixing the coefficients dxy^z- 



D. Dimer coefficients of SU(2)-symmetric models 

At the S'[/(2)-symmetric AF point, the A term in 
the effective Hamiltonian (Q) becomes marginal and in- 
duces logarithmic corrections to several physical quanti- 
ties. Such a logarithmic fashion often makes the accuracy 
of numerical methods decrease. Instead of numerical ap- 
proaches, usinffi-the asymptotic form rf-ithe spinxorre- 
lation functionEZi and OPE technique Jj'EI OrignactJ has 
predicted 
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^ V2 



1/4 



0.2269 



(24) 



at the S'C/(2)-symmetric point. Substituting Eq. ( p4| ) into 
Eq. (pT[), the spin gap in a S'C/(2)-symmetric AF chain 
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FIG. 5: (Color online) Excitation gaps Ai?AS^ =±i (circles) 
and AEas' =o (triangles) of the dimerized XXZ model (H) 
with Az = 0.6. Solid and dashed-dotted curves are fitted 
by the gap formula (|ll[ ) with coefficients outside and inside 
parentheses in Table ^ respectively. The dashed curve repre- 
sents the lightest breather mass which is fixed by combining 
the solid soliton curve and Eq. (lid). 



with dimerization S^y — 5z = S {H 
as 



XXX-5 



AEsu2/J = 1.723(5^/^ 



) is determined 



(25) 



The marginal term however produces a correction to this 
result. It has been shown in Ref. 05^ that the spin gap in 
the model T-p^^^-^ is more nicely fitted with 



AE,^2/J ^ 



1.723,52/3 



1 + 0. 147 In 



0.1616 



1/2 ■ 



(26) 



from the renormalization- group argument. As can be 
seen from Eq. (^6|) , the logarithmic correction is not sig- 
nificantly large for the spin gap. We may therefore apply 



the way based on the sine-Gordon model in Sec. [II C 
even for the present AF Heisenberg model. The resul- 
tant data are listed in the first line of Table fl Evaluated 
coefficients d^y = 0.228 (0.204) and dz = 0.110 (0.097) 
are fairly close to the results of Eq. (p4|). This suggests 
that the effect of the marginal operator on the spin gap 
is really small. We should also note that dxy = "idz is ap- 
proximately realized, which is required from the SU{2) 
symmetry. The numerically calculated spin gap A£'su2j 
Eq. ( |2^), and the curve of the gap formula (|ll|) are shown 
in Fig. B|(a). It is found that even the curve without any 
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FIG. 6: (Color online) Ratio between two numerically eval- 
uated gaps A_Eas^ =o/A_Eas^ =±i (circles) in the dimer- 



ized chain (H) with Sx 



0.3 and dz 



0. Solid curve 



is the soliton- breather mass ratio Ebx/Es in the effective 
sine-Gordon theory (bl). Note that in the ferromagnetic side 
Az < 0, there is no breather and Ebi is replaced with the 
mass gap of soliton-antisoliton scattering states 2Es-, namely, 
EBjEs-f2Es/Es^2. 
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logarithmic correction can fit the numerical data within 
semi-quantitative level. At least, parameters dxy,z at the 
S'C/(2)-symmetric point can be regarded as effective cou- 
pling constants when we naively approximate a dimerized 
Heisenberg chain as a simple sine-Gordon model. 

As discussed in Sec. p Q, logarithmic corrections van- 
ish in the J-J2 model (p^ due to the absence of the 
marginal operator. As expected, Fig. @(b) shows that the 
spin gap Ai?su2 is accurately fitted by the sine-Gordon 
gap formula (O) in the wide range < (5 < 0.3. There- 
fore, the coefficients dxy,z of the J-Ji model (the final line 
of Table E) are highly reliable. Remarkably, the difference 
between the values outside and inside the parentheses is 
much smaller than that of the Heisenberg model (the first 
and last line of Table |). Here, to determine d^y^z of the J- 
J2 model, we have used its spinon velocity v — 1.174Ja, 
which has been evaluated in Ref. pq. 



E. Coefficients of spin operator 



In this subsection, we discuss the spin-operator coef- 
ficient a\ in Eq. (||). Although a\ for the, easy-plane 
XXZ modeliias been evaluated analyticallyclrE3 and nu- 
merically,OL3 those for the S'f7(2)-symmetric Heisen- 
berg chain and the J-Ji model have never been stud- 
ied. The existent data also help us to check the validity 
of our method. From the bosonization formula (0), the 
^-component spin correlation function has the following 
asymptotic form: 



(5|5|,) 



1 A\{-\y'^' 



(27) 



FIG. 7: (Color online) Numerically evaluated gaps 
A^AS" =±1 of dimerized XXZ chains with several values of 
both parameters S^y and Sz at Az = 0.6 and —0.2. Solid 
curves are Eq. (nil) with d^y and dz in Table |. In the fer- 
romagnetic case Az — —0.2, the analytical curve success- 
fully fits the numerical data for a wide weakly-dimerized 
regime 5xy,z ^ 1, while the deviation occurs for the strongly- 
dimerized one. 



in the easy-plane TLL phase. The amplitude Af is re- 
lated to oi as 



Al = al/2. 



(28) 



Lukyanov and his coUaboratorsoO have predicted 
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Jq t \sinh(77t) cosh[(l — 77)<] 
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(29) 



The same amplitude has been calculated by using DMRG 
in Refs. |4|j2|. 

In order to determine ai, we use XXZ models in 
a st agger ed field (fsh. Following the similar way to 
Sec. IIIC| , we can extract the coefficient ai by fitting nu- 
merically evaluated gaps of the model (n3) through the 
sine-Gordon gap formula ([l5|). We numerically estimate 
the gaps at hj J = 0.01, 0.02, •••, 0.09, 0.1, 0.2, and 
0.3 via Aitken-Shanks method. The results are listed in 
column (C) of Table O. Similarly to the case of dimer- 
ization, we adopt spin gaps at relatively large staggered 
fields hsj J = 0.1 and 0.3 to determine the coefficients ai. 



TABLE II: Spin-operator coefficients aiof spin-^ XXZ chain and the J-J2 chain. Values in column (A), (B), and (C) correspond 
the analytical prediction from Refs. [2l|-|23|, the result by DMRG in Refs. [24|j25|, and ours, respectively. 
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soliton gap Es/J 
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0.9 

0.8 
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0.6 
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0.4 

0.3 

0.2 

0.1 



J-J2 model 



0.7049 
0.6069 
0.5464 
0.5008 
0.4629 
0.4297 
0.3994 
0.3712 
0.3443 
0.3183 



0.64 

0.587 

0.54 

0.499 

0.4626 

0.4297 

0.3995 

0.3713 

0.3443 

0.3183 



0.4724 
0.5327 
0.5226 
0.5019 
0.4771 
0.4505 
0.4235 
0.3966 
0.3701 
0.3440 
0.3183 
0.4693 
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(0.4830) 
(0.4808) 
(0.4693) 
(0.4530) 
(0.4338) 
(0.4127) 
(0.3903) 
(0.3670) 
(0.3430) 
(0.3183) 
(0.4668) 
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0.7952 

0.7468 

0.7048 

0.6667 
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0.5970 

0.5641 

0.5319 
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FIG. 8: (Color online) Spin gaps (circles) of (a) the Heisen- 
berg model with d imer ization 5xy = 5z = 5 and (b) the dimer- 
ized J-J2 model (17a). Both solid curves in panels (a) and 
(b) are determined from the gap formula (llll). The dashed 
curve in panel (a) represents Eq. (E6|). 



The value outside (inside) the parentheses in Table O cor- 
responds to fli fixed at hs/ J = 0.1 (0.3). Note that the 
XY model in a staggered field is solvable through Jordan- 
Wigner transformation, and as a result the coefficient ai 
is exactly evaluated as 

ai = I/tt = 0.3183. (30) 

The table clearly shows that the values at hg/ J — 0.1 are 



closer to those of the previous prediction in Refs. 21-25 
We emphasize that our results gradually deviate from 
the analytical prediction from Eq. d29) as the system ap- 
proaches the S't/(2)-symmetric point. The same property 
also appears in the DMRG results in Refs. ^,^. Actu- 
ally, A\ in Eq. (|2^ ) diverges when A^ — >■ 1. However, 
the bosonization formula (0) for spin operators must be 
still used even around A^ = 1. Thus we should real- 
ize that the relation (Eq) is broken and ai remains to 
be finite at the S't/(2)-symmetric point. Figure ^ rep- 
resents the numerically evaluated gaps /S.Ej:^s^^^=±i, and 
three fitting curves fixed by oi (A) and oi (C) outside 
and inside the parentheses in Table Q. Our coefficient oi 
successfully fits the numerical data semi-quantitatively 
in the wide regime 0.01 < hg/ J < 0.3, while the curve 
of fli (A) is valid only in an extremely weak staggered- 
field regime < hs/J < 0.01. This implies that when 
A^ is near unity, the field theory description based on 
Eqs. ( psi) and (|29|) is valid only in a quite narrower re- 
gion for the present staggered-field case compared to the 
case of dimerized spin chain. On the other hand, Fig. y 
also suggests that if we use oi (C) in Table || as the ef- 
fective coefficient of bosonized spin operator instead of 
ai (A) and (B), the XXZ chain in a staggered field dl3) 
may be approximated by a simple sine-Gordon model in 
wide region 0.01 < hs/J < 0.3. 

At the S'f/(2)-symmetric point A^ = 1, a logarithmic 
correction to staggered-field induced gaps is expected to 
appear due to the marginal perturbation. This makes 
it difficult to extract the value ai within the present 
sine-Gordon framework. According to the prediction in 
Ref. |l5| baaed on the asymptotic form of spin correlation 
function^] ai is given by 



1 /7r\i/4 
ai = - - =0.3564 



(31) 



at the S'C/(2)-symmetric point, where ai = bo is imposed. 
The spin gap in AF Heisenberg chains in a staggered field 
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FIG. 9: (Color online) Spin gaps AEASf^^=±i (circles) of 
the XXZ models in a staggered field at A^ = 0.9. Solid 
and dashed-dotted curves are determined from the gap for- 
mula (hSj with ai outside and inside the parentheses in Ta- 
ble |j, respectively. The dashed curve is given by the for- 
mula (|l5|) with ai in Refs. PfJ-feq [i.e., ai determined from 
Eq. mK 



("H'^t'^s with A^ = 1) is thus determined as 



(32) 



A more eorreet gap formula including the logarithmic 
correction has been developed in Refs. p[|lO| as follows: 



A^ 



ASr„. 



,i/J=T85(VJ)2/3[ln(J//.,)]i/6. (33) 



In Fig. ^(a), the numerically evaluated spin gaps, 
Eq. (33), and the fitting curve with ai outside the paren- 
theses in column (C) are drawn. One finds that both 
curves agree well with the numerical data in the weak- 
field regime < hg/J ^ 0.1, while they start to deviate 
from the data in the stronger-field regime. This suggests 
that even at the 5[/(2)-symmetric_point, a simple sine- 
Gordon description for the model (|3|) is applicable in the 
relatively wide region < hs/J < 0.1, if the coefficient 
ai outside the parentheses in column (C) is adopted. 



In the same way as the final paragraph in Sec. [II D, 
we can accurately determine the coefficient ai = bo for 
the J-J2 model since the marginal perturbation vanishes. 
The data are listed in the final line in Table ||. One sees 
from Fig. |l^(b) that the spin gap AEas^^^=±i is fitted 
by the gap formula dlq) quite accurately. In addition, 
the difference between the values outside and inside the 
parentheses is significantly small. 



Coefficients determined from ground-state 
energy 



Instead of the gap formula (|1 1\) , the formula for 
ground-state energy ([l^ ) can also be utilized to deter- 
mine dimer coefficients d^y^z- Let us here define AE'qs = 
Eqs — Ecsi^xy^^z), where E'gs is the ground-state en- 
ergy of the XXZ chain (|]) per site and EQs{Sxy,Sz) is 
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FIG. 10: (Color online) Spin gaps AE^.^^ — n (circles) of (a) 
the Heisenberg and (b) the J-J2 models (]17b|) in a staggered 
field. Solid and dashed curves represent EqTWa) and Eq. (631), 
respectively. We have used ai outside the parentheses in col- 
umn (C) of Table 0. 



that of the bond-alternating XXZ chain (||) . If the dimer- 
ization parameter is small enough li^ajy^^l ^ 1, AEqs is 
expected to agree well with A^qs in Eq. (p^. In this 
case, we can extract the values of dxy,z from the relation 
A£;gs = A^GS- 

To extrapolate the thermodynamic-limit value of 
EGs{Sxy,Sz), we use Aitken-Shanks method for the re- 
sults of finite-size numerical diagonalization, and the 
method works well since the bond-alternating chains are 
gapful. On the other hand, Eqs includes a large finite- 
size correction, as shown in Sec. Ill A. Therefore, instead 
of numerically-ejialuated Eqs , we use its exact value fixed 
by Bethe ansatzo 



E, 



GS 



1 

= - COS 7 - 

4 



1 



• 2 
sm 7 

dA 



cosh(7rA) cosh(27A) — cos 7 



(34) 



where 7 = cos ^ A^. At the limit of 7 — ?> 0, we obtain 
the ground-state energy for the Heisenberg model. 



Eqs 



J 



7-i-O 



= i-ln2. 
4 



(35) 



Black points in Fig. |ll| show Ai^GS determined from 
Eq. (p3) and numerically evaluated EQs{Sxy,Sz) for the 
cases of A^ = 1, 0.9, 0.6 and 0.3. The solid curve in the 
panel (a) of this figure represents the formula (p^) with 
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FIG. 11: (Color online) (a) Ground-state energy difference 
A-Egs for the 5'L''(2)-syininetric case with A^ — 1 and 



Oxy — 0; 

circles) . 
with d^ 



= 6, obtained from numerical diagonalization (black 
Solid and dashed-dotted curves represent Eq. (113) 



!/[^l 



determined from the relation Afcs ~ A_Egs at 



{Sa:y,5^) = (0.05,0) [= (0,0.05)] and with those in Table ^ 
respectively. Dashed curve is Eq. (BGl) including the logarith- 
mic correction. In the panels (b), (c), and (d), black circles 
are A_Bgs for Az — 0.9, 0.6, and 0.3, respectively, under the 
condition 5, = 0. Solid and dashed-dotted curves are respec- 
tively Eq. ( |l2| ) with dxy obtained through A^gs = A_Egs at 
dxy = 0.05 and with that in Table ffl. 



dxy[z] determined from AEqs at {Sxy,Sz) = (0.05, 0)[= 
(0,0.05)]. Solid curves in the panels (b), (c), (d) are also 
the formula (|l2| ) with dxy obtained in the same way. For 
comparison, we also draw dashed-dotted curves of the 
formula (^ with the coefBcients in Table | In the SU{2) 
case of the panel (a), the ground-state energy formula 
with a logarithmic correction 
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0.2728(54/3 
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(36) 



which is predicted in Ref. |15J, is also plotted as a dashed 
curve. As pointed out in Ref. Hq, we find that even the 
curve including the correction deviates from the numer- 
ical data for 6 > 0.1. On top of this isotropic case. 
Fig. |l^ shows that the accuracy of the fitting curves be- 
comes worse as the anisotropy A^ decreases. This is a 
natural result from the fact that the formula (^2|) is bro- 
ken down at the XY point with A^ = and ry = 1/2. 
The deviation between the numerical data and the curve 
also becomes larger for 5 > 0.1 in the easy-plane re- 



gion except for the case around A^ = 0.9. This sharply 
contrasts with the firm correspondence between dimer- 
ization gap and the sine-Gordon gap formula (O) (see, 
e.g. Figs. H"^. We therefore determine the coefficients 
dxy[z] by using the numerical data AE'gs for small dimer- 
ization parameters {Sxy,Sz) = (0.05, 0)[= (0,0.05)] or 
{Sxv,Sz) = (0.1,0)[= (0,0.1)]. They are summarized in 
Table [II. There exists a large difference between d^y^z in 
Tables I and [II , especially, in strongly easy-plane region. 

In the remaining part of this subsection, we discuss the 
reason why AEqs fairly deviate from the analytic pre- 
diction A£gs in contrast to the case of the dimerization 
gap in Sees. IIIB- [II D. Firstly, the sine- Gordon the- 



ory is just a perturbative low-energy effective theory for 
dimerized spin chains, while AEqs would be subject to 
high-energy states as well as low-energy ones. Therefore, 
it is expected that the formula (12) can be applicable 



only in an extremely weak dimerization regime. In fact, 
we find from Fig. O that solid and dashed-dotted curves 
seem to become close to each other in an extremely weak 
dimerization regime Sxy,S < 0.05. Hence, we conclude 
that it is dangerous to apply the sine- Gordon formula of 
the ground-state energy to the original spin chains with 
moderate dimerization. Secondly, the ground-state en- 
ergy difference AEqs is always a convex-downward func- 
tion of Sxy^z in the whole region < A2 < 1. This convex 
property generally makes the accuracy of fitting decrease 
as the case of the dimerization gap in the ferromagnetic 
region A^ < 0. Moreover, as mentioned above, the for- 
mula (n2) becomes invalid in the vicinity of both A^ = 1 
and Az ~ 0. From these arguments, coefficients dxy^z 
and a\ obtained from low-lying excitation gaps are more 
reliable. 



TABLE III: Dimer coefficients dxy^z of spin-i XXZ chain ob- 
tained from ground-state energy difference A_Egs. The data 
outside (inside) the parentheses are fixed by the energy at 



0.05 (0.1). The final line is the result for the J-J2 



model. 



A, 


d 


xy 


d^ 


EGS-EGs{ixy,l^.) 


1 


0.226 


(0.239) 


0.107 (0.113) 


1.148(<5.^d.j, + 5.d.)i-333 


0.9 


0.261 


(0.265) 


0.131 (0.134) 


1.331(<5,„d,„-H5,d,)i*^2 


0.8 


0.275 


(0.274) 


0.143 (0.144) 


\A%^{&xydxy+&AzY-^'''' 


0.7 


0.283 


(0.278) 


0.152 (0.151) 


1.673(5.«d,, + 5,d,)i'503 


0.6 


0.285 


(0.278) 


0.159 (0.156) 


1.924(5,„d,„-H5,d,)i'55o 


0.5 


0.284 


(0.273) 


0.162 (0.158) 


2.2%'d{&xydxy+&AzY-'^ 


0.4 


0.276 


(0.264) 


0.163 (0.157) 


2.827{Sxydxy + Szdz^-''''^ 


0.3 


0.262 


(0.248) 


0.159 (0.132) 


S.766{5xydxy + 5zdz)^'''^° 


0.2 


0.236 


(0.221) 


0.149 (0.140) 


5.705{5xydxy + Szdz)^''"'"^ 


0.1 



J-J2 


0.188 


(0.174) 


0.123 (0.114) 


n.72{Sxydxy + Szdz^-''''' 


0.342 


(0.334) 


0.173 (0.171) 


1.265{5xydxy + SzdzY'''''' 
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FIG. 12: (Color online) Ground-state phase diagram of the 
dimerized spin chains under a magnetic field _ff, Eq. (p7|). 
Each curve represents the phase boundary between the dimer 
and field-induced TLL phases. 



IV. APPLICATIONS 



III 



to some 



In this section, we apply the results of Sec. 
magnetic systems. We demonstrate that several phys- 
ical quantities related to spins or dimerizations can be 
calculated accurately from the data in Tables | and ||. 



A. Dimerized spin chains in a uniform field 



We first consider a spin-^ dimerized XXZ chain in a 
magnetic field. The Hamiltonian is defined as 



U 



S-H 



n 



XXZ-S 



h'^Sj, 



(37) 



with Sxy — 6 and Sz — ^zS. As we have already ex- 
plained, a spin gap opens in the zero-field case. How- 
ever, a magnetic field H > induces the Zeeman split- 
ting, and the gap of the magnon excitation with 5'^ = 1 
(—1) decreases (increases) as AE/^s" =±i T H. When 
H becomes larger than the value of the zero-field spin 
gap, the S*^ = 1 magnon condensation takes place and 
a field-induced TLL phase emerges with an incommen- 
surate Fermi wave number kp = t: — 2t:{S'j). Therefore, 
the curve of the spin gap as a function of dimerization 6 
is directly interpreted as the ground-state phase bound- 
ary of the model y^qj, if the vertical axis (spin gap) is 
replaced with the strength of the magnetic field H. It is 
shown in Fig. |lj. 

The critical point between the dimerized and TLL 
phases can be determined from experiments with varying 
H. Comparing the experimentally obtained critical field 
He and the phase diagram of Fig. n3 in quasi ID dimer- 
ized spin-i compounds, one can evaluate the strength of 
the dimerization S. 



B. Two-leg spin ladder with a four-spin interaction 

We next consider an S'C/(2)-symmetric two-leg spin-i 
AF ladder with a four-spin exchange, whose Hamiltonian 
is given by 

"H *" = 2^ 2^ JSr,j ■ Sr,j+1 + 2_^ J±Sl.j ■ S2,j 
j r=l.2 j 

+ 2^ J4{Sij ■ Sij+i){S2j ■ S2,j+i). (38) 



The symbol r denotes the chain index. Three quantities 
J > 0, J± and Ji respectively stand for the intrachain- 
, interchain- and four-spin coupling constants. There 
are at least two kinds of physical origin of the four-spin 
term J4. The first is that optical phonon modes wr 
a spin-Peierls type coupling can cause a negative J4 
The second is that the higher-order expansion of hop- 
ping terms in half-filled electron ladders with a strong on- 
site Coulomb repulsionoo In fact, the cyclic exchange 
term defined on each plaquette in the ladder contains 
a positive J4 term, which is known to have scaling di- 
mension 1 and be most relevant in all the four-spin cou- 
plings of the cyclic term in the weak rung-coupling regime 

J»|J±|,|J4|. 

The j-jjoadcl (^) has been analyzed by some 
groups .OCJCj There appear four kinds of competing 
phases: the rung-singlet- HaJdane, columnar-dimer, and 
staggered dimer phases.ll3cJ In particular, the ground- 
state phase diagram in the region of J^ > and J4 > 
has been numerically completed in Ref. [44|. 

Here, we show that the data in Tables B and p allow 
us to construct the phase diagram of the model (3^) in 
the weak rung-coupling regime with reasonable accuracy. 
From the bosonization, the low-energy effective Hamilto- 
nian of Eq. (38) reads 



njlad 

"•eff — 






{Sdf 



a 2 2 



1 d^ 
-J±d'^cos{V2TT9- 



Ji^^) cosiVs^cj). 



(39) 



Here we have defined boson fields (f>± — (01 ± (j)2)/V2 
and 6± ~ {di ± 6'2)/\/2, where 0r and Or are dual fields 
of the r-th chain (see Sec. p A[ ). In Eq. (39), we have ex- 
tracted only the most relevant part of the rung couplings. 
The SU(2) symmetry requires the relations v — 7rJa/2, 
K = 1/2, ai = bo = a and d^y = 2dz = 2d. Due to this 
symmetry, three vertex terms in Eq. ( p9| ) have the same 
scaling dimension 1. The {(/)+, 9+) sector is equivalent 
to a sine-Gordon model. A Gaussian-type transition is 
expected at J±d^ — J4{3d)'^ = if other irrelevant pertur- 
bations are negligible. On the other hand, the {(j)-,9-) 
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sector is a self-dual sine-Gordon model,ca which is known 
to yield an Ising-type transition due to the competition 
between cos(-\/87r</'-) and cos(-\/27r^-)- The transition 
occurs as the strength of two coupling constants becomes 
equal, namely, \J^__a? + J4(3d)^|/2 — \J_La?\. Since we 
have already obtained the values of a and d (see Tables || 
and O), we can draw the phase transition curves in the 
J±-Ji space in the weak rung-coupling regime, which are 
shown in Fig. 13, The two transition curves are repre- 
sented as 



/ a 






(40a) 
(40b) 



Each phase is characterized by the locked boson fields 
and their position: In the columnar [staggered] dimer 
phase, (f>+ and (j>- are respectively pinned at \/tt/8 and 
[0 and y^] and (-l)J(5'ij • ^u+i + S2,j ■ S2.J+1) oc 
{sm{V2^cf,+ )cos{V2^^^)) ^ [(~l)^'(5i,j • Si,,+i - 
S2,j ■ S2.J+1} OC (cos(\/27r0+) sin(-\/27r0_)) ^ 0]. In 
the rung-singlet (Haldane) phase, 0- is pinned instead 
of (f>- and ((/>+) = -^/tt/S (0), which corresponds to a 
non-zerO|-ijeMpi"-("odd"-)type nonlocal string order pa- 
rameter .tStZio 

It has been shown in Ref. ^ that Eq. (^) can be 
fermionized. The resulting Hamiltonian consists of three 
copies of massive Majorana fermions and another one 
(For detail, see e.g. Refs. H,0,n6|). The mass of the Ma- 
jorana triplet Mt and that of the remaining one Ms are 
given by 

Mt (xJ^a^ - Ji{Mf, (41a) 



Ms (xSJjLO^ + JiiSdf 



(41b) 



The transition curves in Fig. |3| are identified with Mt — 
and Ms = 0. At Ms = 0, the low-energy physics is 
governed by the gapless singlet fermion which is equiv- 
alent to a critical Ising chain in a transverse field. The 
transition at Ms = therefore belongs to the Ising uni- 
versality class with central charge c = 1/2. On the other 
hand, three copies of massless Majorana fermions, which 
appear at Mt = 0, are equivalent, to an SU{2)2 Wess- 
Zumino-Witten (WZW) theoryErLl with central charge 
c = 3/2. Thus, the transition at Mt = is expected 
to be a c = 3/2 (firi>ifO«i<|;r) type if the marginal current- 



current interactionta'E2IE3 omitted in Eq. (p9) is irrelevant 
(relevant). In Ref. El, the transition has been proved 
to be described by a SU{2)2 WZW theory at least in 
the region of J 3> J±, J4 > 0. This suggests that the 
marginal term is irrelevant there. The Majorana fermion 
with the mass Mt corresponds to a spin-triplet excitation 
(magnon), and another fermion with mass Ms is a spin- 
singlet excitation, which is believed to be continuously 
connected to two-magnon bound state observed in the 
strong rung-coupling regime. 

Finally, we note that in the extremely weak rung- 
coupling limit, the coupling constants of vertex opera- 
tors in Eq. (p9|) would be less valid since coefficients a 



\ 



Ising transition \ 
(c=l/2)\^^;,- 

Haldane 



Stagiered dimer / 




/ c=3/2 or Ist-order 
^ transition 




FIG. 13: (Color online) Ground-state phase diagram of the 
spin ladder (BSl) in the weak rung-coupling regime. There are 
two transition curves, J4 « 2.05Jx and J4 ~ —6.15Jx. The 
former is c = 3/2 or first-order type, while the latter is in the 
Ising universality class with c = 1/2 (see the text). 



and d are determined from gaps induced by relatively 
large staggered field (Jisj J = 0.1 or 0.3) and dimeriza- 
tion (bxy,z — 0.1 or 0.3), respectively. The true transition 
curves might somewhat deviate from our prediction (Hfl) . 
Our result is expected to be more reliable in a moder- 
ate rung-coupling regime. In fact, a numerical study in 
Ref. ^ has shown that the phase boundary is located at 
J4/J± ~ 2 around Jj_/J = 0.25 (see Fig. 6 in Ref. |4|), 
being consistent with Eq. (40a). We stress that our co- 
efficients a and d provides an easy way of estimating 



the phase boundary although it is a rough approxima- 
tion compared with other sophisticated strategies such 
as DMRG and renormalization-group calculations. If we 
replace th e intrachain term in Eq. ( |38|) with two J-J2 
chains 
in Eq. 



), the intrachain marginal interaction omitted 
(391) disappears. In this case, the prediction from 
the effective theory ( |39| ) becomes more reliable even in 
the weak rung-coupling limit J_l/ J, J4/ J — > 0. From the 
data of the J- J2 model in Tables | and O, two transition 
curves in the modified ladder are 



J4 sa0.69J_L, 
J4 « - 2.08 J_L. 



(42a) 
(42b) 



C. Optical response of dimerized spin chains 

Optical responses in Mott insulators including mul- 
tiferroic compounds have been investigated intensively. 
Quite recently, the authors in Ref. ^ have theoretically 
studied the optical conductivity in a ID ionic-Hubbard 
type Mott insulator with Peierls instability, whose strong 
coupling limit is equal to a spin-^ dimerized Heisenberg 
chain, TH?^^^'^ . The results in Ref. E^ would be relevant 
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to, for example, organic Mott insulators such as TTF- 
BA.e3 In this system, the uniform electric polarization P 
along the ID chain is shown to be proportional to the 
dimer operator: 



p^9CiY.^-iysrs,+i, 



(43) 



where g is the coupling constant between the polarization 
and dimer operators. Therefore, P can be bosonized as 



(44) 



P KMg I dxiiin{\/ATT(j){x)) + 



with dxy = 2(i and d 



d. From Eq. (H3) , one can calcu- 
late P and related observables by means of the bosoniza- 
tion for the dimerized spin chain. It has been shown that 
the spin-singlet excitation, i.e., the breather with mass 
Eb^ , is observed as the lowest-frequency sharp peak in 
the optical conductivity measurements. Since the mass 
Eb2 is evaluated from the sine- Gordon theory as 



EbJJ = ViEs/J = 2.924<52/3, 



(45) 



we can extract the value of 6 from the peak position of 
the optical conductivity. The exact expectation value 
of vertex operators in the sine-Gordon model has been 
predicted in Ref. |2|. According to it, the polarization 
density is calculated to be 



{P)/L - {A/3f^'{Esa/vy/^3dg, 



(46) 



with A ~ 3.041 and L being the chain length. This 
provides an experimental way of estimating the coupling 
constant g, which is usually difficult to determine in other 
multiferroic compounds. 



new strategy relying on a solid relationship between the 
lattice models and their low-energy effective theories. 
Our numerical approach is relatively easy compared 
with another method based on DMRG, developed in 
Refs. p4,p5l although the accuracy is expected to be 
better in the latter method. The obtained coefficients 
are summarized in Tables | and || and Fig. |[ In 
addition to these coefficients, we have pointed out a 
dangerous nature of applying the correlation ampli- 
tude ( |29| ) as coefficients of bosonized spin operators 
near the 5f/(2)-symmetric point A^ = 1 in Sec. [HE. 
Furthermore, we have also used the formula for ground- 
state energy of sine-Gordon model to calculate the 
same dimer coefficients in Sec. IIIF| . We conclude that 
the excitation-gap formula (|ll| ) is more suitable than 
the ground-state energy formula (n3) for determining 
coefficients of bosonized operators. 

Physical quantities associated with dimer and spin op- 
erators can be evaluated accurately by utilizing the dimer 
and spin coefficients. As examples, we have determined 
ground-state phase diagrams of dimerized spin chains in 
a uniform field and a two-leg spin ladder with a four- 
spin interaction in Sec. |V|. In addition, we have shown 
how to estimate the electromagnetic coupling constant 
and the strength of the dimerization from the optical ob- 
servables in a ferroelectric dimerized spin chain. These 
applications clearly indicate high potential of the data in 
Tables | and 0. 

An interesting future direction is to apply a similar 
method to other ID systems including fermion and boson 
models. Our method in this paper can be applied to 
lattice systems which have a well-established low-energy 
effective theory, in principle. 



Acknowledgments 



V. CONCLUSIONS 

We have numerically evaluated coefficients of 
bosonized dimer and spin operators in spin-^ XXZ 
model (p and J-J2 model (|l^), by using the correspon- 
dence between the excitation gap of deformed models 
with dimerization (or with staggered Zeeman term) and 
the gap formula for the sine- Gordon theory. This is a 



The authors thank Kiyomi Okamoto, Masaki Os- 
hikawa, and Torn Sakai for useful comments. S.T. and 
M.S. were supported by Grants-in-Aid for JSPS Fellows 
(Grant No. 09J08714) and for Scientific Research from 
MEXT (Grant No. 21740295 and No. 22014016), respec- 
tively. The program package, TITPACK version 2.0, de- 
veloped by Hidetoshi Nishimori, was used in the numer- 
ical diagonalization in Sec. [II. 



^ T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, New York, 2004). 

^ I. Affleck, in Fields, Strings and Critical Phenomena, 
LesHouches, Session XLIX, edited by E. Brezin and J. 
Zinn- Justin (Elsevier, Amsterdam, 1989), p. 564. 

^ V. E. Korepin, N. M. Bogoliubov, and A. G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Functions 
(Cambridge University Press, Cambridge, England, 1993). 

* M. Takahashi, Thermodynamics of One- Dimensional Solv- 
able Models (Cambridge University Press, Cambridge, 



England, 1999). 

A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems (Cambridge 
University Press, Cambridge, England, 1998). 
A. M. Tsvelik, Quantum Field Theory in Condensed Mat- 
ter Physics, 2nd ed. (Cambridge University Press, Cam- 
bridge, England, 2003). 

P. D. Francesco, P. Mathieu, and D. Senechal, Conformal 
Field Theory ( Springer- Verlag, New York, 1997). 
F. C. Alcaraz and A. L. Malvezzi, J. Phys. A, 28, ^521 



15 



(1995). 

M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, ^883 



(1997). 
^° I. Affleck and M. Oshikawa, Phys. Rev. B 60, |l038| ( 1999) 
" F. H. L. Essler and A. M. Tsvehk, Phys. Rev. B 57, 

(1998). 
^^ I. Kuzmenko and F. H. L. Essler, Phys. Rev. B 79, 



10592 



024402 



(2009). 

^^ F. D. M. Haldane, Phys. Rev. B 25, |4925| (1982). 
^'^ T. Papenbrock, T. Barnes, D. J. Dean, M . V. Stoitsov, and 



M. R. Strayer, Phys. Rev. B 68, |024416| (2003). 
^^ E. Orignac, Eur. Phys. J. B 39, 335 (2004). 
^^ D. G. Shelt on. A . A. Nersesyan, and A. M. Tsvehk, Phys. 



Rev. B 53, |852l| (1996). 
^^ E. H. Kim, G. Fat h, J. Solyom, and D. J. Scalapino, Phys. 



14965 



(2000). 



127202 



Rev. B 62, 
^* O. A. Starykh and L. Balents, Phys. Rev. Lett. 93, 

(2004). 
" M. Kohno, O. A. Starykh, and L. Balents, Nat. Phys. 3, 



79C 



(2007). 



077205 



^" OTA. Starykh and L. Balents, Phys. Rev. Lett. 98, 

(2007). 
^^ S. Lukyanov and A. B. Zamolodchikov, Nucl. Phys. B, 493, 

571 (1997). 



^2 S. Lukyanov, Phys. Rev. B 59, |lll63| (1999). 
^^ S. Lukyanov, Nucl. Phys. B 654, 323 (2003). 
^^ T. Hikihara and A. Furusaki, Phys. Rev. B 58, |R583 

(1998). 
^^ T. Hikihara and A. Furusaki, Phys. Rev. B 69 

(2004); cond-mat/0310391 



064427 



^^ F. H. L. Essler and R. M. Konik, arXiv:cond-mat/ |0412421 

^■^ S. Lukyanov, Nucl. Phys. B 522, 533 (1998). 

^* D. C. Cabra, A. Honecker and P. Pujol, Phys. Rev. B 58, 



3241 



B 



(1998). 

. Zamolodchikov, 



Int. J. Mod. Phys. A, 10, 1125 



32 
33 
34 
35 

36 

37 



A 

(1995) 

D. C. Dender, P. R. Hammar, D. H. Rei ch, C. Broholm, 

and G. Aepph, Phys. Rev. Lett. 79, |l75C| (1997). 

K. Okamoto and K. Nomura, Phys. Lett. A 169, 433 

(1992). 

J. L. Cardy, J. Phys. A : Math. Gen., 17, |L38E| (1984). 

J. L. Cardy, Nucl. Phys. B, 270. 18 6 (1986). 

J. L. Cardy, J. Phys. A, 19, |L1093| (1986). 

D. Shanks, J. Math. Phys. (Cambridge, Mass.), 34, 1 

(1955) 



16814 



(1 993) 



K. Nomura, Phys. Rev. B 48, 

I. Affleck, J. Phys. A : Math. Gen. 31, |1573| (1998). 

K. Okamoto and T. Nakamura, J. Phys. A, 30, 

(1997). 



3287 



41 
42 



C. N. Yang and C. P. Yang, Phys. Rev. 150, ^ (1966). 
A. A . Nersesyan and A. M. Tsvehk, Phys. Rev. Lett. 78, 
39391 (1997). 



M. Takahashi, J. Phys. C 10, [l289| (1977). 

A. H. MacDonald, S. M. Girvin, and D. Yoshioka, Phys. 



Rev. B 37, |9753| (1988). 
'^^ A. K . Kolezhuk and H.-J. Mikeska, Phys. Rev. Lett. 80, 
27091 (1998). 

A A Wt spi. ... 



K. Hiju and K. Nomura, Phys. Rev. B 80, |014426| (2009). 

P. Lecheminant, A. O. Gogolin, and A. A. Nersesyan, Nucl. 

Phys. B 639, 502 (2002). 

M. Nakamura, Physica B 329-333, 1000 (2003). 

H. Katsura, M. Sato. T . Furuta, and N. Nagaosa, Phys. 

Rev. Lett. 103, |l77402| (2009). 

F. Kagawa, S. Horiuc hi. M . Tokunaga, J. Fujioka and Y. 

Tokura, Nat. Phys. 6, Km (2010). 



